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To go into the details about the variability of the double quasar SBS 0909+532, 
we designed a monitoring programme with the 2 m Liverpool Robotic Telescope 



in the r Sloan filter, spanning 1.5 years from 2005 January to 2006 June. The r- 
band light curves of the A and B components, several cross-correlation techniques 
and a large number of simulations (synthetic light curves) lead to a robust delay 
AieA = — 49 ± 6 days (la interval) that agrees with our previous results (the B 
component is leading). Once the time delay and the magnitude offset are known, 
the magnitude- and time-shifted light curve of image A is subtracted from the 
light curve of image B. This difference light curve of SBS 0909+532 is consistent 
with zero, so any possible extrinsic signal must be very weak, i.e., the observed 
variability in A and B is basically due to observational noise and intrinsic signal. We 
then make the combined light curve and analyse its statistical properties (structure 
functions). The structure function of the intrinsic luminosity is fitted to predictions 
of simple models of two physical scenarios: accretion disc instabilities and nuclear 
starbursts. Although no simple model is able to accurately reproduce the observed 
trend, symmetric triangular flares in an accretion disc seems to be the best option 
to account for it. 
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1 Introduction 



SBS 0909+532 consists of two co mponents (two-image g r avitationaUy l ensed 



Lubin et al. . 2000; Oscoz et al. 



quasar) sepa r ated by about 1.1" ( Kochanek et al. 1 . 119971 ; iLehar et al. 1 . 12000 



1997a|). In optical frames taken at normal 



seeing conditions with relat ively short exposur e times, the lensing elliptical 



galaxy is undetectable (e.g. lUllan et al. 1 . 120061 ) . Thus, a simple photometric 



model (with only two close point-like sources) is able to describe the whole 
crowded region associated with the quasar c omponents. The first resolved light 



curves of SBS 0909+532 were presented by lUllan et al. I (120061 ). who derived 



accurate fluxes in the R Johnson-Cousins-Bessel filter from a multisite ob- 
serving campaign in 2003. These first i?-band records were used to obtain 
the time delay between both components. From two different techniques and 
1000 repetition s of the expe r iment (synthetic light curves based on the ob- 
served records), Ullan et al. I ( 20061 ) reported a delay At B A ranging from — 41 
to — 56 days (> 90-95% confidence), where the minus sign means that the 
intrinsic signal is observed first in the faintest component (B) and later in the 
brightest component (A). However, the first variability study had some weak 
points. The i?-band light curves did not permit to rule out positive delays 
fairly, and a negative interval [— 90, 0] days was considered in the estimation 
of uncertainties (component B leading component A). Moreover, there was a 
relatively poor overlap between the A and B records, when the A light curve 
was shifted by the best solutions of the time delay and the magnitude offset. 

The light curves of the two components of a double quasar at a given wave- 
length, provide extremely valuable astrophysical information. As the ray paths 
are different for different components, the corresponding traveltimes will not 
agree with each other: it appears a time delay between the observed compo- 
nents, which is directly related to the present expansion rate of the Universe 
(Hubble constan t ) and the mean surfa ce density of the lensing galaxy (e.g. 



Kochanek et al. I . 120041 ; iRefsdal I . Il964j ). To tackle the determination of H 



and < k >, one needs to measure the time delay and the basic parame- 
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ters of the gravitational mirage, i.e., redshifts and image positions, where the 
relevant image positions are the positions with respect to the centre of the 
lensing galaxy (the matter /energy content of the Universe plays also a role). 
Apart from the time delay, a magnitude offset can be also derived from the 
comparison of both brightness records. Thus, the possible extrinsic variability 
(due to intervening objects, e.g., microlenses or dusty clouds in the lensing 
galaxy) may be detected through the difference light curve (DLC): the differ- 
ence between the light curves of the components, w hen the time delay and the 



magn itude offset are taken into account properly (jSchmidt &: Wambsganss 
1998h . While the PLCs of QSO 0957+561 do not show evidence f or ex- 



trinsic variability (] Gil-Me rino et al. I . l200ll ; ISchmidt fc Wambsganss I . 11998 



Wambsganss fc Schmidt I . I1998I ). the DLCs of other lens syste ms seem to indi- 



cate the existence of extrinsic gradients and fluctuations (e.g. iParaficz et al. 
2006h . 



When the difference signal is consistent with zero, the variability observed in 
both components can be attributed to observational noise and intrinsic phe- 
nomena, i.e., physical processes in the source quasar. In this hypothetical case 
ther e is a unique opportunity to study the intrinsic signal of a distant quasar 



e.g. iKawaguchi et al. I . Il998l ). The light curves of the two components can be 



combined to produce one better-sampled record, and this combined ligh t curve 



(CLC) can be used to carry out statistical analyses. IKawaguchi et al. I (|1998 ) 



analy sed the logarithmic slope of the structure function (e.g. lSimonetti et al. 



19851 ) of a CLC of QSO 0957+561, which (slope) is in reasonable agreement 
with the predictions by an accretion disc-instability model. They remarked 
one important caveat about the methodology to discriminate between differ- 
ent scenarios, which is based on the logarithmic slope of the structure function 
of the intrinsic luminosity for nuclear starbursts and accretion disc instabili- 
ties. The anal ysis only incorporated two simple models, i.e., a standard star- 



burst model (jAretxaga et al. 1 . 119971 ) leadi ng to relatively high sl opes and a 



cellular-automaton disc-instability model (jMineshige et al. I . Il994l ) producing 



the smallest slopes, so the conclusions could be biased if the simple models do 
not describe the actual behaviour of the two physical scenarios or the intrinsic 
signal is caused by several independent mechanisms. However, the comparison 
with simple analytical or numerical models represents an important first step 
to reveal the origin of the intrinsic signal. 



In order to get high-quality information about the variability of SBS 0909+532 
at a red wavelength, we designed a monitoring programme with t he world's 



large st fully robotic telescope. The 2 m Liverpool Robotic Telescope (ISteele et al. 



2004 ) at the Roque de los Muchachos Observatory (Canary Islands, Spain) was 
used from 2005 January to 2006 June to obtain nightly frames in the r Sloan 
filter. In this paper we present the robotic programme and the corresponding 
light curves (Sect. 2). Section 3 is devoted to the time delay estimation from 
the records in the r band. Unfortunately, the lens galaxy of SBS 0909+532 
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has a large e f fectiv e radius, with a correspondingly low surface brightness 



( jLehar et al. I . l2000l ). These photometric properties complicate the goal of de- 



termining an accurate galaxy astrometry, and thus, accurate values of H and 
< k >. Even using the more recent optical frames taken with the Hubble 
Space Telescope (HST), the position of the centre of the galaxy has a large 
uncertainty. Consequently, we do not discuss the Hq and < k > values, and 
focus on a careful analysis of the nature of the observed variability. In Sect. 
4, we make the DLC and look for possible extrinsic signal. The CLC and 
its statistical properties (structure function) are analysed and interpreted in 
Sect. 5. Finally, in Sect. 6 we summarise our conclusions and put the results 
in perspective. 



2 Observations and light curves 



The monitoring programme with the 2 m Liverpool Robotic Telescope (LRT) 
began in 2005 January. This robotic project was carried out with the RATCam 
optical CCD camera. The field of view and the pixel scale (binning 2x2) 
were ~ 4.6' x 4.6' and 0.278 arcsec, respectively. Although we concentrate 
here on the observations in the red arm of the optical spectrum, i.e., the 
frames in the r Sloan passband, the gravitationally lensed quasar was also 
monitored in the blue arm (via LRT in the g Sloan passband). A chromatic 
(gVrR) multisite study will be presented further on. Our red subprogramme 
was optimized to get frames all nights when SBS 0909+532 is visible and there 
are no technical/atmospheric problems. Each useful night we usually obtained 
one r-band frame of the lens system (exposure time of 120 s). However, many 
nights in the 2005 October-December period we obtained two 120 s frames, so 
we are able to check formal photometric errors by comparing them with the 
artificial (untrue) intranight variabilities. In fact, the intranight scatters are 
the best non-biased estimators of the typical photometric errors. In Table 1 we 
include details about the whole r-band monitoring from 2005 January to 2006 
June (month, frames and observing nights, sampling rate, and total exposure 
time). Note that this robotic project is quite cheap in observation time, since 
the total science time is only ~ 18 ks = 5 hours. However, this kind of 1.5-year 
monitoring (see Table 1) is organizationally very c omplex or impossibl e with 
a conventionally scheduled and operated telescope ( Steele et al. . 2000l ). 



Basic instrumental reductions are applied to all RATCam frames before the 
data are passed to users. This incorporates bias subtraction, trimming of the 
overscan regions and flat fielding. After the basic pre-processing, cosmic ray 
rejection was applied to the frames. There are als o bad pixel masks, w hich are 



kindly made available by the Angstrom project (IKerins et al. 1 . 120061 ). which 
is another gravitational lensing programme underway on the LRT. Additional 
details on the CCD, the Sloan filters and the data pipeline (LRT automatic 
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Table 1 

Liverpool Robotic Telescope observations of SBS 0909+532 in the r Sloan filter 



Month 


Frames x nights 


Sampling rate (nights/week) 


Total exposure 


2005 January 


1x2 


0.5 


240 


2005 February 











2005 March 


1x12 


3.0 


1440 


2005 April 


1x2 


0.5 


240 


2005 May 


1x7 


1.7 


840 


2005 June 


lxl 


0.2 


120 


2005 July 


occultation 


— 


— 


2005 August 


occultation 


— 


— 


2005 September 


occultation 


— 


— 


2005 October 


lxl + 2x6 


1.7 


1560 


2005 November 


1x5 + 2x11 


4.0 


3240 


2005 December 


1x4 + 2x12 


4.0 


3360 


2006 January 


lxl + 4x1 


0.5 


600 


2006 February 


1x8 + 2x1 


2.2 


1200 


2006 March 


1x16 + 2x1 


4.2 


2160 


2006 April 


lxll + 2x1 


3.0 


1560 


2006 May 


1x13 


3.2 


1560 


2006 June 


1x2 


0.5 


240 



pre-processing) can be found at http://telescope.livjm.ac.uk/ 



A si mple photometric model works well for SBS 0909+532 (see Introduction 
and lUllan et al. 1 . 120061 ). Due to the faintness of the lensing galaxy, the whole 
crowded region can be described through only two close point-like sources 
(two components of the lensed quasar). Thus, we apply a PSF fitting method 
to a representative subset of LRT optical frames verifying some elemental 
conditions (e.g., the telescope pointing was accurate enough so that the lens 
system is included in the field of view, there is no a strongly degraded signal, 
etc) . The subset contains 92 frames, so 60% of the robotic exposures in Table 1 
are used to make the light curves. In Fig. 1 we show subframes corresponding 
to three successful exposures: 2005 March 23 (top panels), 2005 December 4 
(middle panels) and 2006 March 18 (bottom panels). The double quasar (left 
panels) appears as an extended structure similar to two tangent s tars (see the 
right p anels including the "c" field star). Following the notation of lUUan et al. 
( 120061 ). we obtain measurements of yA = mA—mb, Ub = rn>B—rnb, Va = rn a —mb 
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Fig. 1. Liverpool Robotic Telescope subframes of SBS 0909+532 (left panels) and 
"c" field star (right panels) in the r Sloan filter. The top, middle and bottom panels 
correspond to exposures in 2005 March 23 (FWHM = 1.49"), 2005 December 4 
(FWHM = 1.07") and 2006 March 18 (FWHM = 1.15"), respectively. All the 
subframes (panels) have similar size, and the brightness of the "c" star is similar 
to the brightness of the A component. The double quasar (SBS 0909+532) roughly 
looks like two tangent stars. 

and y c = m c — nib, i.e., we compute relative magnitudes using the "b" star 
as the reference object (the lens system is inside the triangle denned by the 



a-c" stars, e.g. iKochanek et al. 1 . 119971 ). Both "a" and "b" are non-variable 



field stars, and we expect a constant behaviour of y a . To check the trend of 
y a , we focus on the 69 relative fluxes corresponding to the 2005/2006 season: 
continuous monitoring from 2005 October to 2006 June. Using the formal 
photometric errors, the trend is not consistent with a constant because the 
reduced chi-square value is large (x 2 = 6.27). The formal uncertainties do not 
seem to include all the sources of error, and we need a non-biased estimator of 
uncertainties. This is not a problem in our project, since we have several nights 
with two exposures at different times. Taking into account the deviations in 
y a for nights with two data (intranight deviations: 5k = y a (tk + At) — y a {tk), 
At < 8 hours, k — 1, N), it is easy to obtain a standard intranight deviation 



Ca = \Y^k=i3l/(N — 1) = 6.8 mmag. Assuming this standard deviation as a 
typical error (i.e., the intranight variability is not real, but the tool to estimate 
the true photometric uncertainty) and re-doing the fit to a constant, we derive 
a very reasonable \ 2 value of 1.12 (see the filled squares and the discontinuous 
line after day 3600 in Fig. 2). Therefore, the m a -m;, record in the 2005/2006 
season is clearly consistent with a constant behaviour and there is a fair way 
to infer non-biased uncertainties of the quasar fluxes y& and y% (typical errors 
from standard intranight deviations). 

From intranight variabilities in the 2005/2006 season, we derive self-consistent 
uncertainties: <J a ~ 7 mmag < a c ~ 10.5 mmag < <ja ~ 14 mmag < <jb ~ 
18 mmag. We remark that the "c" star and the A component have a similar 
brightness, but A is placed in a crowed region and reasonably <ja is larger 
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3400 3600 3800 

JD-2450000 



Fig. 2. Liverpool Robotic Telescope light curves of SBS 0909+532 in the r Sloan 
filter. The filled circles are the fluxes jja, shifted by + 0.50 mag, and the open circles 
are the fluxes ys (see main text). We also incorporate relative fluxes of the "a" star 
for comparison purposes. The filled squares represent the y a data, shifted by + 
3.13 mag, and the discontinuous line is the fit to the fluxes corresponding to the 
2005/2006 season (after day 3600). While the stellar fluxes are narrowly distributed 
around a constant flux, the A fluxes trace gradients and have a larger error, and 
the B fluxes have the largest uncertainty and scatter. 




3400 3600 3800 

JD-2450000 

Fig. 3. Grouped light curves of SBS 0909+532. We grouped the individual quasar 
fluxes (see Fig. 2) within 3-day intervals. The filled circles are the fluxes da, shifted 
by — 49 days and + 0.65 mag (see section 4), and the open circles are the fluxes 

VB- 

than a c . While the stellar relative fluxes are inferred with true uncertainties 
less than or equal to 10 mmag, the quasar relative fluxes have true errors in 
the 10-20 mmag range. Once non-biased errors in y^ and ys are determined 
from intranight variations, we group the quasar fluxes each night including 
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more than one exposure. In the 2005/2006 season, the quasar light curves are 
characterized by a mean sampling rate of two points every week. The LRT was 
offline in 2005 September (for two weeks) for engineering work which included 
realuminisation of the primary mirror. Thus, our whole original light curves 
(incorporating data before 2005 September) are affected by this maintenance 
work, and we must correct the realuminisation offsets, i.e., the instrumental 
offsets between 2004/2005 and 2005/2006 seasons. There is no problem with 
y a , since we can easily estimate the difference between the average fluxes before 
and after day 3600. However, the quasar components are variable objects, and 
we need independent brightness records to correct the offsets in both y^ and 
yB- A simultaneous monitoring in the R band at Mt. Maidanak, Uzbekistan, 
is used to derive the realuminisation biases in the LRT records. Maidanak (1.5 
m AZT-22 Telescope) quasar fluxes are compared with LRT brightnesses at 
separations less than or equal to 3 days, which leads to instrumental offsets 
of about 70 and 10 mmag in y^ and yB, respectively. As the R — I and B — R 
colours of the "b" star are close to the colours of the B component (similar 
spectra), it is expected a very small correction in the differential curve ys- 
For this record ys, we just obtain an offset (slightly less than 10 mmag) 
that agrees with theoretical predictions. The whole final brightness records in 
the r band are depicted in Fig. 2. In this figure, we plot the A light curve 
(filled circles), shifted by + 0.50 mag, and the B light curve (open circles). To 
compare quasar fluxes against stellar brightnesses, relative fluxes of the "a" 
star (filled squares), shifted by + 3.13 mag, are also shown in Fig. 2. The LRT 
brightnesses in the r-band (78 data of each component) have errors above 10 
mmag, which could complicate some analyses. Thus, to improve the situation 
we also group fluxes within 3-day intervals. The new data set only contains 41 
fluxes of each component, but both mean errors are reduced to the 10-mmag 
level: < oa >~ 10 mmag and < ob >~ 13 mmag. In Fig. 3 we show the 
new grouped light curves, i.e., the A curve (filled circles) and the B curve 
(open circles). Now, to easily compare both trends, the A record is properly 
shifted in time and magnitude (— 49 days and + 0.65 mag; see sections 3-4). 
The individual calibrated quasar fluxes are presented in Table 2. We use the 
SDSS flux of the "b" star in the r band (m& = 14.87 mag) to get apparent 
magnitudes of both component 



3 Confirmation of the time delay 



If we concentrate our attention in Fig. 3, the A and B light curves show a 
decline larger than 50 mmag as well as a 50-mmag event around day 3800 
(A is shifted in time and magnitude). Thus, although there is an important 
gap in the LRT monitoring (due to occultation of the lens system in 2005 



The SDSS Web site is http://www.sdss.org/ 
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July-September), these two features are promising tools to confirm our pre- 
vious time delay estimation (see Introduction). To calculate the time delay 
between both components of SBS 0909+532, we use different cross-correlation 
techniques. First, we focus on the grouped light curves and the discrete cross 
ation function (DCF). The DCF was introduced by lEdelson fc Krolik 



corre 



2002 



(|1988) and has been extensively used in delay studies (e.g 



Gil-Merino et al. 



Ofek fc Maoz~ll2003l : IOscoz et al. I . ll997bh . Here, the DCF of SBS 0909+532 



is evaluated every 1 day in the region from — 200 to + 100 days, so we analyse 
a wide range of lags including both positive and negative values. The DCF 
is binned in 2a day intervals centered at the lags. To work with a reasonable 
time resolution, we only take into account a values less than or equal to 10 
days, i.e., bins with width < 20 days. For a < 4 days, the DCF is very noisy, 
whereas for a = 8-10 days, the main peak of the DCF is significantly reduced 
with respect to the expected value of 1. The more interesting results are de- 
rived from a = 5-7 days, and we choose the intermediate bin (a = 6 days) 
as the most suitable one. For a = 6 days, there is a maximum at — 50 days 
(DCF max ~ 0.9). Apart from the main peak around the maximum (delay- 
peak), there are other secondary peaks at negative and positive lags. These 
secondary structures have an amplitude of about 0.5, i.e., they are clearly 
smaller than the main feature. No correlation (DCF ~ 0) or anticorrelation 
(DCF < 0) is also found at the edges of the time lag-interval. For other a 
values, we also find maxima at or around — 50 days. 



As it was discussed by iLehar et al. I (jl992|) , when the main fluctuations in the 



light curves have an intrinsic origin, the irregular delay-peak of the AB cross- 
correlation function should be closely traced by the symmetrical central peak 
(around a lag equal to zero) of the AA (or BB) autocorrelation function. More- 
over, other features of the cross-correlation function around lags Ti, T2,... will 
be closely reproduced in the autocorrelation function around lags T x — AtBA, 
T 2 — At B A,---, respectively (with relation to A, we assume that B is delayed 
in At ba)- Therefore, if the shifted discrete autocorrelation function (DAF) 
is matched to the discrete cross-correlation function (DCF), one derives the 
time delay in a self-consistent way (note that the lag corresponding to the 
maximum of the DCF is only a rough estimation of the delay). This self- 
consistent methodology is called the S 2 technique, and it has been successfully 
applied to some golden data sets (see next paragrap h) for QSO 0957+561 (e.g. 



Goicoechea et al. . 1998; Serra-Ricart et al. . 199 



apn 

3). 



Before applying the 5 2 technique, we need to make a golden data set. A golden 
data set contains a variable and free from long gaps light curve A in a period 
[ti, tf], and a variable and free from long gaps light curve B in a period [ij+At, 
tf + At], where At ~ AtBA- As we know a rough initial estimation of the SBS 
0909+532 delay (through the maximum of the DCF, see here above), we can 
use this value (At = — 50 days) to correct the unsuitable 50-day edges. It is 
also clear that the records in Fig. 3 are variable. Therefore, the existence of long 
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Fig. 4. Comparison between the DCF (filled circles) and the DAF shifted by — 
50 days (open circles). The DAF is the average of the AA and BB autocorrelation 
functions. We use a quasi-golden data set (see main text) and a = 6 days, where a 
is the semiwidth of the bins. 



gaps is the only difficulty to make a golden data set for SB S 0909+532. Can we 
relax t he condition on the absence of gaps?. Simulations by lSerra-Ricart et al. 
(119991 ) indicated that the 5 2 method works even in the presence of relatively 
long gaps. Moreover, the main problem with gaps is the possible bias between 
the differences Va~ <Va> and ys— < y B >, which are the key pieces in the 
cross-correlation. In the periods of interest, the two light curves are sampled 
in different ways, so the direct averages of the fluxes < yA >d and < ys >d 
could lead to a bias between the differences y^— < yA >d and jjb— < Vb >d- 
Thus we can obtain a quasi-golden data set provided that < yA > and < y B > 
are carefully determined. Instead of the direct averages of the grouped fluxes 
yA and y B in the periods of interest, the mean values < yA > and < ys > are 
inferred from a method that corrects for sampling bias. 



Once we have a quasi-golden data set, a comparison between the correspond- 
ing DCF (filled circles) and DAF (open circles) is plotted in Fig. 4. The DAF 
is the average of the AA and BB autocorrelation functions, and it is shifted 
by — 50 days, i.e., the rough estimation of the delay. There is no important 
distortions in the features of the DCF as compared with the features in the 
DAF, which strengths the realibility of a delay close to — 50 days. Using a 
= 6 days, possible values of the time delay (9) versus the associated 5 2 {6) 
values normalised by its minimum value 5 2 (8q), are also plotted in Fig. 5. The 



S 2 (9 ) function is define d in E q. (1) and Eq. (7) of iGoicoechea et al. I ( 19981 ) 
and Serra-Ricart et al~l ( 19991 ). respectively. In Fig. 5, a relatively narrow peak 
centered on — 49 days (best value of the delay) is derived. To obtain un c ertain - 



ties, we follow an approach similar to that described by lUllan et al. I (120061 ) . 



We make 1000 repetitions of the experiment, apply the 5 2 minimization (a 
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-100 -80 -60 -40 -20 



(days) 

Fig. 5. Normalised 5 2 function from LRT data in the r band. The S 2 function is 
made from the DCF and DAF in Fig. 4. 




200 



Fig. 6. Correlation coefficient associated with the MCF method. We use the indi- 
vidual fluxes in Fig. 2 and 10-day bins in the component A. 

6 days) to each quasi-golden synthetic data set, and thus obtain 1000 best 
values of the delay. Through the distribution of delays, our la measurement 
is AtBA = — 49 ± 7 days (69.8% confidence interval). This S 2 result confirms 
the previous delay determination from Calar Alto and Maidanak frames in 
2003. Moreover, the new measurement is very robust, since it is not a pre- 
conditioned estimation (a wide range from — 200 to + 100 days is tested) and 
there is a significant overlap between the A and B records, when the A light 
curve is shifted by the best solution of the time delay (see Fig. 3). 



We a l so apply the modifie d cross-correlation function (MCF) technique (IBeskin fc Oknyanskij 



19951 ; lOknyanskij 1 . 119971 ). The MCF combines properties of both sta ndard 



cross-correlation functio ns: the interpolated c ross-c orrelation function by lGaskell fc Sparke 



( 119861 ) and the DCF bv lEdelson fc Krolik I (I198SI ). This time we analyse the 
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individual (non-grouped) fluxes in Fig. 2, and after doing some tests, we use 
10-day bins in the component A. When the MCF is applied to our data in the 
lag interval [— 200, + 200] days, the maximum correlation coefficient (R ma x — 
0.878) corresponds to a lag T max = — 48 days. The correlation coefficient at 
different lags (days) is shown in Fig. 6. We also check if the removal of some 
fluxes from the B record can influence the estimation of (Rmax, T max ). The 
stability test is done in a blind way, i.e., we do not choose by eye the points 
to be dropped, but a systematic procedure is used to select those points. In 
a first step, the first 3 points are removed from the B light curve, then the 
second, third and fourth points are removed from, etc. From this 3-point sys- 
tematic cleaning, almost all the R max values are close to 0.9 and ~ 85% of 
the iterations lead to T max = — 49 days, in good agreement with the analysis 
of the whole data set. In order to estimate delay errors, we carry out simu- 
lations taking into account all the properties of the observed records : kind of 



variability, sampling and photometric errors (IKoptelova et al. 1 . 120061 ) . About 



1000 synthetic data sets leading to R ma x > 0.8 are considered here. Our la 
measurement is AtsA = — 49 ± 5 days (70.6% confidence interval), which is 
practically identical to the la result from the 5 2 method. We consider that 
the two techniques (5 2 and MCF) have similar quality. Thus, there is no a 
fair way to choose either 5 or 7 days as error, and we adopt — 49 ± 6 days as 
the final la estimation, where ± 1 day is the uncertainty in the error. 



4 Difference light curve 



From now on we take the LRT grouped fluxes in Fig. 3 as basic tools for 
discussing variability properties in the r filter. First, considering the time delay 
in section 3, we infer the time delay-shifted light curve of A (y^ s )- Here, y\ 
is the da record shifted by AtsA = — 49 days in the horizontal direction. Then 
a mean offset < yB — y T A S > is computed in a direct way. This magnitude offset 
is close to + 0.65 mag. Second, the magnitude- and time-shifted light curve of 
image A (y^ TS , where y^ TS is the y\ s curve shifted by the magnitude offset 
in the vertical direction) is subtracted from the light curve of image B (ys)- 
Both y^ TS and ys are plotted in Fig. 3. We thus obtain the difference light 
curve (DLC) ys — yA TS ■ To compute the mean offset as well as the DLC, the 
dates in the time shifted curves y\ s and y^ TS are taken as reference epochs. 
The y^ s and yjf fluxes are then compared to the averaged values of ys in 
bins with semiwidth a centred on the reference dates. 

The DLC of SBS 0909+532 in the years 2005-2006 appears in Fig. 7. We use 
two different a values: 6 days (filled circles) and 15 days (open circles), and 
both trends are consistent with each other. The differences (in r-magnitudes) 
do not exceed the ± 0.05 mag thresholds (discontinuous lines). Moreover, there 
is no evidence in favor of the existence of events or gradients. The difference sig- 
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Fig. 7. Difference light curve of images B and A (shifted in time and magnitude) 
of the lensed quasar SBS 0909+532. To obtain the differences in r -magnitudes, we 
use bins in B with semiwidths a = 6 days (filled circles) and a = 15 days (open 
circles). The ± 0.05 mag thresholds are depicted by discontinuous lines. 

nal is in apparent agreement with zero, i.e., Fig. 7 shows a noisy relationship ys 
— Ua TS - From a quantitative point of view, we can also estimate the goodness 
of rep resenting the data with y B — y^ TS = (e.g. ISchmidt fc Wambsganss 



19981 ). The reduced chi-square values are 0.92 (a = 6 days) and 0.87 (a = 
15 days), corroborating the good agreement between the light curves of the 
two images. Therefore, the observed variability in A and B is basically due 
to observational noise and intrinsic signal, although we cannot rule out the 
existence of a very weak extrinsic signal whose amplitude must be well below 
the noise level in the DLC. Our constraints on the possible microlensing (ex- 
trinsic) variability can be used to obtain information on the granularity of the 
matter in the lensing g alaxy and the size of the source (Gil-Merino et al., in 
preparation). Although iMediavilla et al. I (120051 ) suggested the existence of an 
achromatic microlensing magnification of image B, a homogeneous microlens- 
ing pattern will not produce microlensing variations but a microlensing offset 
as part of the magnitude offset. Only an inhomogeneous microlensing pattern 
can be detected through the DLC. 



5 Combined light curve: origin of the intrinsic signal 



The DLC presented in Sect. 4 is in clear agreement with the absence or 
an extremely low (undetectable) level of extrinsic signal. Therefore, in this 
section, we make the combined light curve (CLC) and interprete it as due to 
the observational noise and intrinsic phenomena. The combined photometry 
consists of both light curves ?/^f T5 and ys, and this CLC is plotted in Fig. 3 
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(filled and open circles). The data points in Fig. 3 show the global record y qso 
(in r-magnitudes) in a 500-day interval, with the data covering a period P ~ 
450 days (there is an unavoidable gap of about 50 days that is related to the 
annual occultation of the quasar). 

The CLC is generated by intrinsic signal (s) and observational noise (n), so 
Uqso = s + n. This simple relationship between the combined fluxes, the under- 
laying intrinsic signal and the nois e permits to estimate the first -order struc- 



ture function of s = y qso — n (e.g. iGil-Merino et al. 1 . 1200 ll ; ISimonetti et al 



19851 ). A structure function analysis is a method of quantifying typical flux 
variabilities at different lags. The structure function SF(s) at lag At is given 
by 



SF(s) = (l/2iV) - y qso ,) 2 - o\ - a) 



where the sum only includes the pairs verifying that tj — U ~ At (the 
number of such pairs is N). We take a normalization factor equal to 1/2, and 
thus, the asymptotic b ehaviour on long timescal es is just the signal variance 



a 2 instead of 2a 2 (e.g. ICollier fc Peterson 1 . 120011 ). Fig. 8 shows the structure 
function of SBS 0909+532 obtained from the CLC and Eq. (1), using indepen- 
dent 10-day bins (filled circles). To test the robustness of the derived SF(s), 
we also obtain extreme structure functions, i.e., using the extreme values of 
the delay range, At B A = — 43 days and At B A = — 55 days, to make the 
corresponding CLCs. These two extreme trends (continuous lines in Fig. 8) 
are included in the error bars of the SF(s), so the typical structure function 
(from AtsA — — 49 days) is a reliable tool. The discontinuous line in Fig. 8 
corresponds to a 10 mmag threshold, and at lags At < 60 days, the typical 
fluctuations are below this 10 mmag level. This result explains why the flux 
ratio at the same time of observation basically coincides with the flux ratio 
corrected by the tim e delay, at least in the red arm of the optical spectrum 



flUllan et al. 1 . 120061 ). 



To discuss the origin of the intri nsic signal, we focus on the s tructure function 
of the intrinsic luminosity L (e.g. lCid Fernandes et al. 1 . 120001 ) . This SF(L) can 
be directly compared to the predictions of differen t physical scenarios, e.g ., 



nuclear starbursts and accretion disc instabilities (IKawaguchi et al. I . Il998l ). 
We remark that the observed flux at A is emitted at a shorter wavelength 
A/(l + z s ), and we must consider the emission of near UV light (~ 2600 
A), because the observations were made in the r Sloan band (A ~ 6200 A). 
On the other hand, s = m qso — = — 2.5 \og(F qso / F b ), where m and F 
denote r-band magnitudes and monochromatic fluxes, respectively. We can 
also use the cosmological law F qso = eL qso /[A'nD 2 L {l + z s )], with L qso , e and 
Dl being respectively monochromatic luminosity of the quasar, extinction- 
magnification factor and luminosity distance. From some rearrangement, it is 
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Fig. 8. First-order structure function of the underlaying intrinsic signal s in the 
r band. The SF(s) from the CLC in Fig. 3 (filled circles) is compared with two 
extreme structure functions associated with CLCs for the extreme values of the 
delay range (continuous lines) . A discontinuous straight line represents the 10 mmag 
threshold. 
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Fig. 9. Structure function of the intrinsic luminosity L. Taking into account the 
asymptotic behaviour (dotted line), the total monitoring period and the number of 
pairs in the bins, we distinguish between valid/reliable results (filled circles) and the 
rest of data points (open circles). The SF(L) at Ar > 63 days (rest-frame lags on 
the right of the arrow) is probably biased, and thus, the results at long rest-frame 
lags should be considered as unreliable data. While the solid line is the best fit from 
symmetric triangular flares in an accretion disc, the dashed line traces the best fit 
from standard nuclear starbursts (see main text). 
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inferred the relationship kL qso = 10 _a4s , k = e/[47rD|(l + z s )F b ]. We initially 
take a set of units so that k = 1 and L = L qso = lCT a4s , so the structure 
function SF(L) at rest-frame lag At = At/(l + z s ) can be estimated through 
the averaged sum 

SF(L) = (1/2N) ^[(1(T - 4 ^ - 10- - 4 *'".*) 2 - a 2 - a 2 }, (2) 



where cr = 0.921 x 10 0A yi so <j and the sum includes N pairs verifying Tj — Tj ~ 
Ar. As the [SF(L)] 1 / 2 values from Eq. (2) are ~ 10 , we decide to use more 
appropriate units: k = 1/1000. In these new units, the structure function of 
the intrinsic luminosity is drawn in Fig. 9 (filled and open circles). In this 
figure, the dotted straight line represents the asymptotic behaviour ox, i.e., 
[SF(L)] 1 / 2 — > <7x at long lags. From now on, we only consider the data points 
not exceeding the flat asymptotic behaviour (filled circles). There is an arrow 
to make mark on the last reliable point. This last valid point corresponds to 
At ~ 150 days ~ P/3. Moreover, at At < 150 days, the number of pairs in 
each bin is iV ~ 100-200. Longer lags cannot be seriously considered in the 
analysis, since the associated information is probably biased and could not 
describe the typical behaviour of the signal. 

The reliable trend in Fig. 9 (filled points) should be related to the variability 
scenario and may unveil the origin of the intrinsic fluctuations. First, we use 
three simple analytical models of accretion disc instabilities. In these three 
Poissonian models, the luminosity is due to the superposition of a variable 
component and a constant background (the nonflaring part of the accretion 
disc). The variable component is made by the superpositions of flares at ran- 
dom times, which are characterized by a typical timescale T. Square (SQ), 
exponentially decaying (ED) and symmetric triangular (ST) flares are taken 



into account (see Appendix B of ICid Fernandes et al. 1 . 120001 ) . When the ob 



served SF(L) is fitted to the analytical intrinsic structure functions, the SQ 
and ED models do not work well and lead to x 2 ~ 4-5. However, the ST model 
leads to a more reasonable (but not sufficiently good) \ 2 value of 2 (solid line 
in Fig. 9). For the best ST model, T = 70 days and the lifetime is r = 3T/2 ~ 
100 days. There are several physical timescales that might be associated with 
accretion disc instabilities. The main dynamical timescales are the free-fall 
time Tf f ~ (FP/GM) 1 ' 2 and the orbital time r orb ~ 2tt(R 3 /GM)^ 2 ~ 6r f/ . 
The thermal timescale is given by r t h ~ r ff/li i- e - ; the time for vertical dif- 
fusion of heat. Apart from the black hole mass M and the emission radius R, 
the thermal timescale also depends on the disc viscosity parameter 7 (this last 
parameter is usually named a, but we rename it 7 to avoid confusion). The 
viscous time r visc ~ (R/h) 2 r th (h is the disc thickness) is also relevant. This is 
linke d to the radial i nflow of a gas element, and represents t he longest timescale 



(e.g. ICzernv 1 . 12004 ; iKrolik I . Il999l . and references therein) . ICollier fc Peterson 



( 120011 ) estimated the values of Tff, r orb and r th for different values of M at 
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two rest-frame wavelengths: 1400 and 5000 A (see Fig. 5 in that work). To 
do the est imations, they assum e d a st andard geometrically thin and optically 
thick disc IShakura Sz Sunyaev I (119731 ). At ~ 2600 A (see above), the lifetime 
of ST flares is consistent with the orbital time for ~ 10 9 M Q black holes and 
the thermal timescale for ~ 10 7 M P) black holes. 



Many previous studies focused on the logarithmic slope ((3) of the growing part 
of the square root of structure function, which (slope) is directly related to 



the p hysical mechanism responsible for the variability (e.g. iKawaguchi et al. 



19981). For example, for SQ flares of duration T, [SF(L)] 1 ' 2 oc (At) 1 / 2 at 



lags Ar < T (SF(L) is flat at Ar > T). This leads to a constant slope (3 
= 0.5. For ED and ST flares, (3 reaches its maximum value at the shortest 
lags. If the ED flares are characterized by a semi-lifetime T, then (3 ~ 0.5 at 
Ar < T/2. On the other hand, for ST flares with rise time T, [SF(L)] 1 ^ 2 is 
approximately proportional to Ar at Ar < T/2. Thus, whereas (3 should be 
less than or equal to 0.5 for the SQ and ED models, the ST model produces a 
steep slope (3 ~ 1 at relatively short lags. What about the observed structure 
function?. The observed slope is (3 ~ 1 using both the subset of data at Ar < 
35 days (the first six points in Fig. 9) and the global data set until Ar = 
63 days. It is clear that the observed slope of about 1 favours the ST model 
of accretion disc instabilities. However, even ST flares with equal rise and 
decay times T = 70 days (best solution) have not the ability to reproduce 
the observed slope at lags exceeding 50 days (see Fig. 9). In fact this is the 
reason to obtain a rough agreement (x 2 = 2) instead of an accurate fit leading 



to x ~ 1- The cellular-automaton model (IMineshige et al. I . Il994j ) is another 



framework to describe accretion disc instabilities. In spite of its popularity, 
we cannot consider that physical mechanism as an alterna tive to explain the 



obser vations, since it generates a small slope: (3 = 0.4-0.5 (IKawaguchi et al. 



19981 ) . 



With respect to the nuclear starbursts, we use the standard model by lAretxaga et al. 



(119971 ). This simple model is characterized by a timescale r sg , which is the time 
when the supernova remnant reaches the maximum of its radiative phase. 
The observed SF(L) c annot be reproduced by the relationship in Eq. (17) of 



Aretxaga et al. I (119971 ). Although the best timescale has an acceptable value 
of r sg = 60 days (considering that the source is a very bright quasar), the 
observations-model comparison leads to a poor fit with x 2 = 4 (dashed line 
in Fig. 9). Standard nuclear starbursts are able to induce a steep slope at lags 
shorter than r sg . However, the flat regime is only reached at lag s significantly 



longer than r sg (lAretxaga et al. 1 . 119971 : IKawaguchi et al. I . Il998l ) . This mecha- 



nism is thus unable to reproduce a steep slope from short lags to the flattening 
lags, and consequently it fails to trace the observed behaviour. 
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6 Conclusions and discussion 



The 2 m Liverpool Robotic Telescope ( Steele I . l200ll ; ISteele et al. 1 . 120041 ) at 
the Roque de los Muchachos Observatory (Canary Islands, Spain) is ideally 
suited to monitorize gravitationally lensed quasars (GLQs) and derive light 
curves of their components. In this paper we present the first GLQ monitoring 
campaign using the Liverpool Telescope. We have observed the double quasar 
SBS 0909+532 in 2005-2006, taking nightly frames in the r Sloan filter. Our 
main results and conclusions are: 



The r-band frames and PSF photometry lead to variable light curves of 
the two quasar components A and B. Both brightness records show a 
decline larger than 50 mmag followed by a 50-mmag event. These fea- 
tures are use d to confirm the re cently reported time delay between the 



components (jUUan et al. 1 . 120061 ). From different cross-correlation tech- 



niques and a large number of repetitions of the experiment (synthetic 
light curves based on the observed records), we infer a delay AtsA = — 
49 ± 6 days (la interval), which agrees with our previous results. The 
new delay determination is robust, since a wide range of possible delays 
is tested and the two light curves overlap over a long time interval, when 
the A record is shifted by the best solution of the time delay. 
(2) To obtain the difference light curve of SBS 0909+532, the magnitude- 
and time-shifted light curve of image A is subtracted from the light curve 
of image B. There is no evidence in favor of the existence of extrinsic vari- 
ability (e.g., fluctuations or gradients caused by microlensing in the lens- 
ing galax y halo), since t h e diffe rence curve is consistent with zero. Very 



recently, iParaficz et al. I (120061 ) derived difference curves of five GLQs 
and found that two out o f the five GLQs have significant extrinsic sig- 
nal. IParaficz et al. I (120061 ) also presented one difference curve that clearly 
agrees with zero (HE 2149-2745) as well as two doubtful cases. Here we 
report on another GLQ with flat difference curve, which joins the fam- 



ily of systems without important extrinsic variations (e.g. iBurud et al. 



20021 ; iGil-Merino et al. I . l200ll ; ISchmidt fc Wambseanss I . Il998h . 



(3) Taking into account the absence or the extremely low (undetectable) level 
of extrinsic signal, we make the combined light curve and interprete it 
as due to the observational noise and intrinsic variations. The combined 
photometry consists of both A and B records, where the A light curve is 
shifted by the best solutions of the time delay and the magnitude offset. In 
order to study the growth of intrinsic variability with rest-frame lag, we 
explicitly obtain the structure function of the intrinsic luminosity. This 
is then fitted to predictions of simple models of accretion disc instabili- 
ties and nuclear starbursts. With respect to the dis c-instability scenario. 



we con centrate on the phenomenological mod els bvlCid Fernandes et al. 



( ]2000l ). since the cellular-automaton model (IMineshige et al. I . 1 19941 ) is 
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related to a relatively sm all logarithmic slope of the structure function 



(IKawaguchi et al. I . Il998l ). This small slope ((3 = 0.4-0.5) is not consis- 
tent with the structure function derived from o bservations (0 ~ 1). W e 
also consider the standard nuclear starbursts by Aretxaga et al. I ( 1997 ). 



Symmetric triangular flares in an accretion disc lead to the best (but not 
accurate enough) fit, whereas standard SN explosions cannot produce a 
large value of (3 from short to long lags, and thus, cannot account for the 
observed variability. 

We obtain a rough agreement between the observations of SBS 0909+532 
and the production of 100-day symmetric flares in the accretion disc of the 
distant source quasar. Considering typical supermassive black holes, the life- 
time of the flares seems to agree with the dynamical a nd thermal times at 



the e mission radius corresponding to a standard disc (ICollier fc Peterson 



20011 ) . This suggest the existence of disc instabilities having the local (emis- 
sion ring) dynamical-thermal timescale. Local magnetorotational instabilities 



e.g. iBalbus fc Hawley I . Il99ll ) would be possibly able to cause variability over 



this timescale. However, it is difficult to identify precisely the physical ori- 
gin of the fluctuations, and we cannot rule out other processes. For example, 
the hottest temperature kT at the innermost radius of a standard gas disc is 
clearly smaller than X- r ay energies. Therefore, a s SBS 0909+532 is a bright 



X-ray source (j Chart as 1 . 12000| ; iPage et al. 1 . 120041 ). inverse Comptonization in 



a hot inner corona may explain its X-ray emission. Instabilities in the corona 



which is presumably unstable (e.g. IShakura fc Sunyaev I . Il976l ). could lead to 



variable X-ray emission as well as variable X-ray irradiation of the gas disc. 
After X-ray r eprocessing; in t he disc, variability at op tical /UV wavelengths is 



expected (e.g. ICzerny 1 . 120041 ) . iManmoto et al. I (119961 ) also proposed an unsta- 



ble advection-dominated disc that might generate X-ray flares and subsequent 
optical/UV events. In both physical schemes, the shape and timescale of the 
optical/UV variations are determined by the shape and timescale of the X-ray 
fluctuations. Are these reasonable pic tures?. X-ray symmetric flares were pre- 
viously found by lNegoro et al. I (119941 ). but on a stellar scale. Cygnus X-l is an 
accreting black hole system (M ~ 10 M^) that produc es symmetric flares last- 



ing about 1 second (see Fig. 1 of lNegoro et al. I . I1994J ). The Cygnus X-l vari 



ability was reproduced through theoretical calculations by IManmoto et al. 



(119961 ). so it may be associated with advection-dominated disc instabilities. If 
the same mechanism is responsible for the X-ray variability of both stellar and 
supermassive black hole systems, then the stellar timescale must be rescaled 
by 10 6 -10 8 (assuming typical supermassive black holes with M ~ 10 7 -10 9 
M Q ) to obtain r ~ 10-1000 days symmetric flares on a quasar scale. For M ~ 
10 8 M , these hypothetical X-ray events would be good tracers of our near 
UV fluctuations. Thus, surprisingly, we find that advection-dominated disc 
instabilities could be able to account for variations on both stellar and quasar 
scales. An unstable hot corona is also a plausible source of X-ray flares, but 
a discussion on the involved timescale and (a) symmetry is out of the scope of 
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this paper. 



To date only one more GLQ with flat difference curve was analysed in some 
detail to determine the origin of its intrinsic variability. Fro m light curves 



of Q SO 0957+561 in two different optical bands (gr filters), iKundic et al 



( 119971 ) inferred variability functions [SF{m)\ 1 ^ 2 and their corresponding log- 
arithmic slope j3 gr ~ 0.4. Here, m = s + n are quasar magnitudes, i.e., the 
observational noise was not subt racted from the p h otome tric measurements. 
Using only data in the g band, 
shallow slope of [SF(m)] 1 / 2 : (3 g r 



Kawaguchi et al. I (119981) a l so det ermined a 



0.3-0.4. iKawaguchi et al. I (119981 ) presented 



a growing trend at lags At < 500 days (At < 200 days), but they did not 
report on the asymptotic behaviour (\[2a m for the Kawaguchi et al.'s normal- 
ization) or the reliability of the results at long lags, so we cannot properly 
discuss the involved variability timescale. On the other hand, the observed 
rise in the variability is consiste nt with a cellular-automaton disc-instability 
model (IKawaguchi et al. I . Il998l ). However, in spite of the small errors in the 
light curves (er ~ 10 mmag), the noise in the structure function (\plo ~ 15 
mmag) does have a significant effect on the measured variations at the shortest 
lags and the measured sope. Thus, the slope of [SF(s)] 1 ^ 2 should be steeper 
than the slope of [S'.F(m)] 1 / 2 . If we roughly construct a noise-less function 
[SF(s)} 1/2 in the g band, the new slope is ~ 0.7. Moreo ver, this new (and 



proba bly true) slope agrees with the slope in the R band (jGil-Merino et al. 



20011 ) . Therefore, there are clear evidences that both QSO 0957+561 and SBS 
0909+532 are z = 1.4 bright quasars with (3 ~ 0.7-1. 



We may also try to establish a relation between these GLQ results and vari- 
ability studies of non-lensed quasars. While some works focused on ensem- 
ble structure fu nctions of large s amples (hundreds, thousands or even more) 



of qu asars (e.g. Ide Vries et al. I . 120051 ; lHawkins I . 120021 ; IVanden Berk et al. 



2000 



2004) , other works cocentrated o n well-sampled individual quasars (e.g. lCid Fernandes et al. 



Collier fc Peterson 1 . 120011 ). The first ones did not include details on in- 



dividual objects. Moreover, additional problems usually complicate a direct 
comparison with current GLQ data. For example, de Vries et al. I ( 2005 ) stated 
that their study is insensitive (due to the measurement noise) to lags shorter 
than t ~ 1 year . Thus, we cannot put the behaviour of GLQs in perspective. 



Hawkins I (120021 ) mostly used relatively faint objects showing coherent fluctu- 
ations of about 1 mag over a few years (see Fig. 4 in that work). However, 
both GLQs are relatively bright sources displaying coherent flux variations of 
about 0.1 mag over a few or several months. On timescales of a few years, the 
GLQs have a maximum scatter (difference between maximum and minimum 
flux) below 0.25 mag. Moreover, instead of a noise- less structure fu nction 
having a good time resolution at lags At < 1-2 years, lHawkins I (120021 ) stud- 
ied a long term ensemble [SF(m)] 1 ^ 2 with time resolutio n of about one year. 



The m ost relevant ensemble variability was reported by IVanden Berk et al. 
(120041 ) . who showed an [SF(s)] 1 ^ 2 vs. At relationship having reasonable time 
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coverage and resolution. This ensemble variability for quasars with very dif- 
ferent luminosities (absolute magnitudes) and redshifts led to a slope (3 ~ 0.3 
at At < 500 days. It is evident that the GLQ slope (3 > 0.5 disagrees with 
the ensemble slope ((3 < 0.5), which could indicate the existence of different 
populations of intrinsically variable quasars (perhaps depending on luminosity 
and redshift) o r even the existe nce of microlensing effects on a large amount 



sources, 



of quasa rs (e.g. [H awkins , 2002 ). Wit h respect to the well-sampled individual 



Cid Fernandes et al. I (I2000h analysed the light curve s of 42 z < 0.4 



quasa rs that were monitored by the Wise Observatory group (IGiveon et al 
19991 ). Assuming the presence of square flares (/3 = 0.5), they derived life- 



times ranging from 124 days to a few yea rs, so the shortest on e s are close 



to the lifetime of the SBS 0909+532 flares. ICid Fernandes et al. I (12000( 1 also 



remarked the insensitivity of their fits to the flare shape (with the excep- 
tion of exponentially decaying flares, since these asymmetric flares led to poor 
fits). These results (relatively long lifetimes and poor sensitivity to the flare 
shape) could be due, at least in part, to the procedure for fitting observations, 
which involved two free parameters, i.e., the flare lifetime and the variance. 
We, however, compute directly the variance from the light curve and then 
fit the lifetime. Our variance represents the typical variability in a period of 
about 200 days (rest-frame time), and taking into account the absence of 
important gradie nts over 1.5 years (from a comparison with previous Calar 
Alto light curves lUllan et al. Ll2006h. this seem s a good tracer of the variance 
over longer periods. ICollier fc Peterson I (120011 ) studied 13 local active galac- 
tic nuclei with very good time coverage and resolution. The slopes of their 
optical/UV structure functions varied between (3 ~ 0.3 and (3 ~ 0.8 (similar 
to the GLQ slope). The flare lifetimes (using symmetric triangular flares and 
certain lag intervals) are r ~ 5-94 days, with the longest ones very close to 
the lifetime of the SBS 0909+532 events. In order to check the results based 



on a monitoring period of 200-300 days, ICollier fc Peterson I (120011 ) presented 
the variability structure for ~ 7 years of monitoring of NGC 5548. Using the 
short monitoring period, they found r ~ 40-60 days. However, from the ad- 
ditional analysis of the longer monitoring period, they inferred richer results. 
The asymptotic behaviour is reached at lags of ~ 200 days, and the 50-day 
timescale seems to represent the lifetime of the shortest flares. 



To accurately reproduce the observed variability of SBS 0909+532, new (more 
complex) models seem to be required. Moreover, a longer monitoring period 
would lead to a very precise description of the structure at At < 60 days 
and the appearance of reliable features at lags exceeding At ~ 60-70 days, 
which would significantly improve the observational constraints. We also note 
that the current optical data of SBS 0909+532 are ins ufficient to ca r ry out 



more ambitious analyses, e.g., a test for non-linearity. lUttley et al. I (120051 ) 



showed that a linear variability-flux relation in the light curve of an object is 
an indicator of non-linearity. This hypothetical linear relation would suggest 
that the variability process is multiplicative, or in other words, it is the result of 
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several random subprocesses which multiply together. A process produced by 
multiplication of many independent processes has a lognormal distribution, 
so the hypothetical light curve would have a logno rmal probability density 
function (PDF) (for details, e.g. lUttley et al. 1 . 120051 ). However, the observed 
fluctuations in the flux of SBS 0909+532 (< 10%) and the number of available 
data points (clearly insufficient to make highly populated bins in flux) do 
not permit to reliably obtain the variability-flux relation or the PDF of SBS 
0909+532. Finally, we remark the unique advantages of using GLQs as a tool 
to study the origin of the intrinsic signal of quasars, because there is no way 
to disentangle intrinsic from extrinsic signal in a non-lensed quasar. New 
long term monitoring programmes of GLQs with flat difference curve will 
permit to fairly discuss the structure of the intrinsic variability on short and 
long timescales, and to compare between this behaviour and the structure of 
the fluctuations in both GLQs with non-flat difference curve and non-lensed 
quasars. 
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Table 2: LRT fluxes of SBS 0909+532A,B in the r Sloan 
filter 



JD-2450000 


rriA (mag) 


<j(7TIa) (mag) 


ma (mag) 


aim p>) (mat 


3446.388 


16.337 


0.014 


16.985 


0.018 


3453.428 


16.336 


0.014 


16.984 


0.018 


3455.432 


16.349 


0.014 


16.977 


0.018 


3458.426 


16.342 


0.014 


16.981 


0.018 


3487.371 


16.333 


0.014 


16.964 


0.018 


3489.371 


16.321 


0.014 


16.970 


0.018 


3503.383 


16.312 


0.014 


16.978 


0.018 


3509.389 


16.325 


0.014 


16.965 


0.018 


3510.384 


16.303 


0.014 


16.989 


0.018 


3522.463 


16.318 


0.014 


16.961 


0.018 


3525.425 


16.317 


0.014 


16.969 


0.018 


3646.719 


16.380 


0.014 


17.060 


0.018 


3649.711 


16.383 


0.010 


17.067 


0.013 


3656.713 


16.395 


0.010 


17.047 


0.013 


3662.709 


16.386 


0.010 


17.051 


0.013 


3667.711 


16.382 


0.014 


17.074 


0.018 


3676.662 


16.395 


0.010 


17.025 


0.013 


3678.723 


16.390 


0.014 


17.032 


0.018 


3680.637 


16.392 


0.014 


17.026 


0.018 


3681.645 


16.391 


0.014 


17.044 


0.018 


3682.660 


16.383 


0.014 


17.047 


0.018 


3684.670 


16.392 


0.010 


17.052 


0.013 


3693.615 


16.390 


0.010 


17.042 


0.013 


3699.588 


16.396 


0.010 


17.051 


0.013 


3700.680 


16.398 


0.014 


17.006 


0.018 


3707.590 


16.399 


0.014 


17.007 


0.018 



Table 2: LRT fluxes of SBS 0909+532A,B in the r Sloan filter (continuation) 



JD— 2450000 rriA (mag) o^m^) (mag) m# (mag) cr(ms) (mag) 



3708.576 


16.390 


0.010 


17.048 


0.013 


3709.699 


16.394 


0.010 


17.030 


0.013 


3711.699 


16.389 


0.014 


17.050 


0.018 


3712.621 


16.389 


0.014 


17.044 


0.018 


3714.590 


16.397 


0.010 


17.048 


0.013 


3718.555 


16.396 


0.010 


17.037 


0.013 


3719.615 


16.391 


0.010 


17.039 


0.013 


3728.629 


16.402 


0.014 


17.037 


0.018 


3729.504 


16.386 


0.014 


17.056 


0.018 


3731.544 


16.394 


0.008 


17.031 


0.010 


3766.609 


16.405 


0.014 


17.030 


0.018 


3783.531 


16.403 


0.014 


17.014 


0.018 


3784.539 


16.399 


0.014 


17.018 


0.018 


3787.457 


16.401 


0.014 


17.025 


0.018 


3789.449 


16.396 


0.014 


17.027 


0.018 


3790.422 


16.400 


0.014 


17.024 


0.018 


3793.410 


16.392 


0.014 


17.015 


0.018 


3798.445 


16.406 


0.014 


16.987 


0.018 


3802.406 


16.390 


0.014 


17.009 


0.018 


3803.395 


16.391 


0.014 


17.002 


0.018 


3804.418 


16.391 


0.014 


17.018 


0.018 


oon r on 1 

3805.391 


16.402 


0.014 


1 t r\r\ a 

17.004 


0.018 


3806.383 


16.387 


0.014 


17.009 


0.018 


3813.395 


16.382 


0.014 


17.019 


0.018 


3814.426 


16.382 


0.014 


17.014 


0.018 


3815.383 


16.381 


0.014 


17.009 


0.018 


3817.438 


16.388 


0.014 


17.034 


0.018 



Table 2: LRT fluxes of SBS 0909+532A,B in the r Sloan filter (continuation) 



JD— 2450000 rriA (mag) o^m^) (mag) m# (mag) cr(ms) (mag) 



3818.426 


16.373 


0.014 


17.032 


0.018 


3820.379 


16.384 


0.014 


17.024 


0.018 


3821.492 


16.376 


0.014 


17.039 


0.018 


3822.492 


16.390 


0.014 


17.028 


0.018 


3826.492 


16.379 


0.014 


17.012 


0.018 


3827.520 


16.380 


0.014 


17.038 


0.018 


3828.492 


16.364 


0.014 


17.032 


0.018 


3830.473 


16.384 


0.014 


17.018 


0.018 


3833.500 


16.368 


0.014 


17.036 


0.018 


3841.465 


16.354 


0.014 


17.063 


0.018 


3843.430 


16.368 


0.014 


17.058 


0.018 


3845.430 


16.367 


0.014 


17.030 


0.018 


3854.402 


16.382 


0.014 


17.028 


0.018 


3855.488 


16.370 


0.014 


17.034 


0.018 


3856.402 


16.384 


0.014 


17.029 


0.018 


3860.391 


16.397 


0.014 


17.026 


0.018 


3862.395 


16.386 


0.014 


17.035 


0.018 


3864.453 


16.378 


0.014 


17.037 


0.018 


3867.438 


16.380 


0.014 


17.041 


0.018 


3870.434 


16.379 


0.014 


17.042 


0.018 


3875.418 


16.372 


0.014 


17.057 


0.018 


3876.477 


16.362 


0.014 


17.036 


0.018 


3879.430 


16.378 


0.014 


17.030 


0.018 


3880.414 


16.361 


0.014 


17.039 


0.018 


3888.418 


16.369 


0.014 


17.025 


0.018 



27 



